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Abstract 



In this paper we prove that the shape optimization problem 



mm 



{A fe (ft) : n C R d , fl open, P(O) = 1, < +00} , 



has a solution for any k € N and dimension d. Moreover, every solution is a bounded 
connected open set with boundary which is C 1,Q outside a closed set of Hausdorff 
dimension d — 8. Our results are more general and apply to spectral functionals of 
the form /(Afc 1 (f2), . . . , Afc p (fi)), for increasing functions / satisfying some suitable bi- 
Lipschitz type condition. 



1 Introduction 

A shape optimization problem is a variational problem of the type 



where A is an admissible family (or admissible set) of domains in M. d and F is a given 
cost functional. For a detailed introduction to this type of problems we refer to the books 
IU, [23 E6] and the papers [HE3]. 

A particularly interesting and widely studied class of shape optimization problem is the 
one in which the admissible set A is composed of open domains with perimeter or measure 
constraint and the cost functional F depends on the solution of some partial differential 
equation on Q. For example, F(p,) = Afe(O), where A&(J1) is the kih eigenvalue of the 
Dirichlet Laplacian, i.e. the kth smallest positive real number such that the equation 



mm{F(Q) : O E A} 



(1.1) 



has a non-trivial solution. Another typical example is 




where p G L 2 (M. d ) and w is the solution of 



Aw = p, 



w e H^(Q). 



Unfortunately, there is no available technique to directly prove the existence of a solution 
of (1.1) in the family of open sets. The standard approach is to extend the cost functional 



to a wider class of domains and then prove that the minimizers are open sets, hopefully 
with a smooth boundary. 

The existence in the case when the domains are measurable sets is well-known if the 
admissible family of domains consists of sets which satisfy some geometric constraint which 
assures a priori some compactness. Under some monotonicity and semi-continuity assump- 
tions on the cost functional F a general theorem of Buttazzo and Dal Maso states that 



there are solutions of (1.1), in each of the cases 

A = {f2 C D : ft Lebesgue measurable, |f2| < c} , (1-2) 

(see [El 16]) and 

A = {O C D ; ft Lebesgue measurable, P{fl) < c, < +00} , (1.3) 

(see [9]). Here |0| is the Lebesgue measure of Q, -P(fi) the De Giorgi perimeter of O (see 
|23| I30|) and the design region D is a bounded open set in M. d . Clearly, in each of these cases 
the functional F has to be appropriately extended on the class of Lebesgue measurable sets 
in D (see Section [2] for more details). 

The case D = M rf is more involved and only recently was proved (see [7] and |31j ) that 



when F(fl) = A&(f2), for some k G N, there exists a solution of the problem (1.1), where 



the admissible set A is given by (1.2) 



Once the existence of a solution of (1.1) is established in the class of measurable sets 
it is quite natural to ask whether these optimal sets are open (hence being also a solution 
of the same problem in the more "natural" class of open sets) and, in case the answer is 
affirmative, which is the regularity of the boundary of these optimal sets. 

The study of the regularity of the optimal set in the case of a measure constraint strongly 
depends on the nature of the cost functional. Even the openness is a difficult question which 
is known to have a positive answer in the case of the Dirichlet energy F(£l) = E p (Q) (see 
[I]) and some spectral functionals F(Q) = /(Ai(O), . . . , A/c(f2)) (see [S], [13] and (T2"]). 



In the cases F(fl) = E p (Q), F(fl) = Xi(fl) the problem (1.1) with admissible set 



(1.2) can be written in terms of a single function in Hq(D), i.e. the corresponding shape 



optimization problems are equivalent to 

mini ^ J \Vw\ 2 dx-J pwdx: w 6 H%(D), \{w ^ 0}\ < cX , (1.4) 



^2 

and 

min 



f f \Vu\ 2 dx ~] 

{ J D J dx : u G ff oP). I{« + 0}| < cj , (1.5) 

respectively. For these functionals one can apply the techniques introduced in [2] to study 
the regularity of the optimal domains (see [1] and [5], for two different arguments based 
on this idea). When the cost functional F(Q) depends on the spectrum of the Dirichlet 
Laplacian the regularity of the optimal sets is still not known in the case of a measure 
constraint. The main difficulty in this case is due to the fact that the variational formulation 
of Afc(fi) does not concern functions but /c-dimensional spaces of functions, which makes 
the analysis quite difficult. More precisely, we have 

f |Vn| 2 dx 

X k (n) = min max ^ — ^- — , (1.6) 
KcH^(n) ueK J n u 2 dx 

where the minimum is over all /c-dimensional subspaces K of Hq{VL). 
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In this paper we are interested in the study of existence and regularity for the problem 
(1.1) under the perimeter constraint (1.3) when the design region is D = M. d . Our main 



result is the following: 

Theorem 1.1. The shape optimization problem 

{A fe (0) : Jlcl^ open, P(O) = 1, |0| < oo} , (1.7) 



mm 



where Afc (Q) is defined in (1.6) through the classical Sobolev space Hq(Q) on the open set 
J7, has a solution. Moreover, any optimal set Q is bounded and connected. The boundary 
dQ, is C 1,a , for every a G (0, 1), outside a closed set of Hausdorff dimension at most d — 8. 



This result is a consequence of the more general Theorem 1 1 . 2 1 which applies to spectral 
functionals of the form 

F(n) = f(\ kl (n),...,x kp (n)), 

where / : W — > R is an increasing locally Lipschitz function satisfying some local bi- 
Lipschitz-type assumption. More precisely, we consider / such that: 

(/l) f(x) — > +oo as \x\ —7- +oo; 

(/2) / is locally Lipschitz continuous; 

(/3) / is increasing, i.e. for any x = (x±, . . . , x p ) G W and y = {y\, . . . , y p ) G W such that 
x > y, i.e. satisfying xj > yj, for every j = l,...,p, we have f(x) > f(y). More 
precisely we assume that for every compact set K cM. d \ {0}, there exists a constant 
a > such that for any x,y G K, x > y, 

f(x) - f(y) > a\x - y\. 

For example, any polynomial of X kl , ■ ■ ■ ,\k v with positive coefficients satisfies (/l), (/2) 
and (/3). 

Theorem 1.2. Suppose that f : W — > R satisfies the assumptions (/l), (/2) and (/3). 
Then the shape optimization problem 



mm 



{/(Aj fcl (n),...,Afc p (n)) : QcR d , n open, P(Q) = 1, |0|<oo}, (U 



has a solution. Moreover, any optimal set Q is bounded and connected and its boundary d£l 
is C 1,a , for every a G (0, 1), outside a closed set of Hausdorff dimension at most d — 8. 

Proof. The existence of an optimal set is first proved in the class of measurable sets in 



Theorem 3.9 The existence of a solution $1 of (1.8) is proved in Theorem 4.8 The 



boundedness of f2 is due to the facts that 17 is an energy subsolution (see Definition 



3.4 



and Proposition 3.5) and that every energy subsolution is a bounded set (Lemma 3.7). 
The problem of regularity of the optimal set is treated in Section [5] More precisely, by 
Proposition |3.5| and Lemma |4.3[ we have that £1 is an energy subsolution and a perimeter 
supersolution. In Theorem |5.6[ we prove that any f2 satisfying those two conditions has 
C 1,a boundary, for every a G (0, 1), outside a closed set of Hausdorff dimension at most 



d — 8. The connectedness of the optimal set follows by Proposition 5.8 □ 



Remark 1.3. The regularity of the free boundary proved in Theorems |1.1| and 1.2 is not 



in general optimal. Indeed, it was shown in [9j that the solution Q of (1.7) for k = 2 



has smooth boundary. The proof is based on a perturbation technique and the fact that 
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A2(0) > Ai (SI) and can be applied for every k G N under the assumption that the optimal 
set is such that Afc(Q) > Afe_i(0), see Remark 5.7 On the other hand it is expected (due 
to some numerical computations) that the optimal set Q for A3 in M? is a ball and, in 
particular, As(f2) = A2(f2). 



Remark 1.4. The bound on the diameter of the optimal set 0, solution of (1.8), depends 
on the function / and on A& (O). In dimension two, this bound is trivially uniform, since 
it depends only on the perimeter constraint. This fact (together with he convexity of the 
minimizers) was used in [11] to study the asymptotic behaviour of the optimal sets f2&, 
solutions of (1.7). More precisely, it was proved that the sequence C R 2 converges in 
the Hausdorff distance (up to translations) to the ball of unit perimeter. The analogous 
result in higher dimensions is not known yet. 



2 Preliminaries 

In this section we introduce the notions and results that we will need in the rest of the 
paper. As we saw in the introduction, we will have to solve partial differential equations 
on domains which are not open sets. For this purpose we extend the notion of a Sobolev 
space to any measurable set Q C introducing the Sobolev-like spaces Hq(Q) defined as 

flo( n ) = { u G ff X ( Rd ) : u = a.e. on O c } , (2.1) 

where the term almost everywhere (shortly a.e.) refers to the Lebesgue measure | • | on M d . 

Remark 2.1. Note that even for open sets f2 the above definition differs from the classical 
one, in which Hq(Q) is the closure of the smooth functions with compact support C^°(Q) 
with respect to the norm 

IMIffi : = IMl£a + Il v vlli2. 

To see that, one may take for example f2 to be the unit ball minus a hyperplane passing 
through the origin. Nevertheless, we have equality Hq(Q) = Hq(Q), if O is a bounded open 
set with Lipschitz boundary or, more generally, an open set satisfying a uniform exterior 



density estimate (Proposition 4.7). 



Remark 2.2. In the case of measure constraint another definition of a Sobolev space is used. 
Indeed, for any measurable set Q C W 1 one may consider 

ffl(O) = {«£ fr 1 (R d ) : cap({n ^ 0} n n c ) = 0} , 

where the capacity cap(E) of a generic set E C M d is defined as 

cap(E') = inf {||u||#i : u G H 1 ^), u = lona neighbourhood of e\ . (2.2) 

In Theorem 3.3.42] it was proved that the above definition coincides with the classical 
one on the open sets of R d . There is a close relation between the Sobolev spaces as defined 
in (2.2) and Sobolev-like spaces from (2.1). It fact, for every measurable set ft C M d , we 
have the inclusion Hq(Q) C Hq(Q,). Moreover, since Hq(Q.) is separable, it is not hard to 
check that there is a measurable set [/Cl li such that U C fl a.e. and Hq(U) = Hq(Q) 
(take, for example, U to be the union of the supports of Sobolev functions, which form a 



dense subset of Hq(£1)). The reason to work with definition (2.1) instead of (2.2) is that 



we do not know the relation (if any) between the perimeter of f2 and the perimeter of U. 
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For any ft C M. d of finite measure and any / £ L 2 (R d ), we define Rn(f) S ^o(^) as the 
weak solution, in #q (fi), of the equation 

-Au = f, u€Hq(Q), (2.3) 
or, equivalently, the unique minimizer in Hq(Q) of the convex functional 

Jf(u) = - I |V«| 2 dx — I fudx. 
Applying the Sobolev inequality in ]R d and using Rn{f) as a test function in (|2.3[), we have 



\\Rn(f)\\h < M 2/d \\Rn(f)\\ 2 ^ < c d \n\ 2/d \MW))\\h < c d \n\ 2/d \\Rn(f)\\ L 4f\\L^ 

L d-2 

and so, R n : L 2 (M d ) -»■ L 2 (M d ) is a bounded linear symmetric operator such that 
\\Rn\\c(L*m) < C d \n\ 2/d , Rn(L 2 (R d )) C ^(fl). 



Moreover, i?^ is compact and positive and so, it has a discrete spectrum cr(Rn) contained 
in M + , which we write as 

< ••• < A fc (0) < < Ai(O). 

We define A/%(f2) as the inverse of Afc(O). Moreover, we have the variational characterization, 
analogous to the one in (1.6): 

~ , , To |Vn| 2 dx 

X k (n) = min max Jn ' , (2.4) 

KcHjtn) ueK Jn u dx 

where the minimum is over all fc-dimensional linear subspaces K of Hq(U). 

For any measurable set f2 of finite Lebesgue measure, we denote with wn the weak 
solution of the equation 

-A«j n = l, w Q eH*(n), (2-5) 
i.e. the rojj is the unique minimizer in Hq(Q) of the functional 

J(w) = - I \Vw\ 2 dx- I wdx. (2.6) 
2 Jn Jn 

We define the energy functional E(£l) as 

E(£l) = min J(w) = J{wn) = -\ I' wndx, (2.7) 
w&Hl{n) ^ Jn 

where to obtain the last equality we have used the definition of J and the relation 

\Vwn\ 2 dx = / wn dx, 



which follows after testing (2.5) with wq. 
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Remark 2.3. For the solution of (2.5) we have the estimates 

IKIIz* < c d \n\ 1+2/d , livwnllia < c d \n\ 3/4+2/d - (2.8) 

Moreover, by [34] , wn G L°° (R d ) and 

\\wn\\ L oo < C d \Sl\ 2/d , (2-9) 

where C d is a constant depending only on the dimension. We also note that in the framework 
of Sobolevdike spaces, the weak maximum principle still holds (wjj < wq, whenever U C 
0), while the analogous of the strong maximum principle is the following equality (see 



Proposition 4.7) 



H^{w a > 0}) = H%({w n > 0}) = H^(Q), (2.10) 



i.e. one may take U = {w^ > 0} in Remark 2.2 



The relation between the operator Rq and the solution , is explained in the following 
two propositions 

Proposition 2.4. Suppose that fl n C M d is a sequence of measurable sets of uniformly 
bounded Lebesgue measure and suppose that WQ n converges to some w G H 1 (M. d ) strongly in 
L 2 (JR. d ), as n — > oo. Then, there is a compact self-adjoint operator R on L 2 (U. d ) such that 
i?n„ converges to R in the strong operator topology in C(L 2 (R d )) and R < Rq in sense of 
operators, i.e. 



R{f)fdx < [ Rn(f)fdx, for every f G IS 



(R d ), 



where we set £1 = {w > 0}. 

Proof. See [U Proposition 3.3]. □ 

Proposition 2.5. Suppose that Q\ and Q2 ire two sets of finite measure in M d such that 
Oi C r^2 ■ Then, there are a positive constant C and a real number 9 G (0, 1) ; depending on 
the measure of 0,2 and the dimension d, such that 

ll-R^i - Rn 2 \\c(L 2 (R d )) < CWw^ - wn 2 \f L2 ^ 

Proof. See [H Lemma 3.6]. □ 



3 Existence of generalized solutions 

In this section, we prove that for every k G N and every c > the problem 

: n C R d , n measurable, P(O) = c, \U\ < ooj , (3.1) 



mm - 



has a solution (see Theorem 3.8). More generally, in Theorem 3.9 we prove that there is a 
solution of 

min{/ (A fel (Q), . . ., A fep (J2)) : QcR d , n measurable, P(fl) = c, \U\ < 00} , (3.2) 

where / is an increasing function satisfying (/l), (/2), (/3). We will use a combination of a 
concentration-compactness principle and an induction argument, as in [7\. In order to deal 
with the dichotomy case, we will show, on each step of the induction, that the optimal sets 
are bounded. The following theorem is a straightforward adaptation of [BJ Theorem 2.2] 
and already appeared in [Sj . We report the detailed proof in the Appendix for the sake of 
completeness. 
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Theorem 3.1. Suppose that Q n C is a sequence of measurable sets of finite measure 
and uniformly bounded perimeter. Then, up to a subsequence, one of the following three 
situations occurs: 

(i) Compactness: There is a sequence y n G M. d , a measurable set ft C M d and a bounded 
self-adjoint operator R : L 2 (R d ) -> L 2 (R d ), such that l Wn+ a n -> In in ^(W 1 ), 
R Vn+ n n ^Rin C(L 2 (R d )), and R < R n . 

(ii) Dichotomy: There are sequences of measurable sets A n C M d and B n C M. d such that: 

(a) A n UB n C n n ; 

(b) d(A n , B n ) — > oo, as n — >• oo; 

(c) liminf^oo \ A n \ > and liminf^oo |Z? n | > 0; 

(d) limsup^ (P(A n ) + P{B n ) - P(Q n )) < 0. 

(e) lim n ^oo \\RA n uB n - RnJ\c(L 2 {R d )) = °- 

(Hi) Vanishing: For every e > and every R > there exists N G N such that 

sup sup \Q n n _Br(x)| < e. 



Moreover, \\Rn n \\c(L 2 (R d )) ~> «« n ^ oo. 
Remark 3.2. Notice that if, in the previous theorem, we assume 



supAi(O ra ) < +oo, (3.3) 

neN 

then the vanishing cannot occur. Indeed, due to the equality 

ll^n„IU(i2 (Rd)) = ~ , 

Al{\l n ) 

the sequence of resolvents does not converge to zero in norm. 

Remark 3.3. If the compactness occurs for the sequence Q n C M. d , then we have 

P(0) < liminf P(tt n ), A fe (0) < liminf A fe (Qj, 

n— >oo n— >oo 

for every fc£N. If the dichotomy occurs, then we have 

liminfP(^ n U5 n ) < liminf P(Sl n ), liminf X k (A n U B n ) < liminf X k (Q n ), 

n— >oo n— >-oo n— >oo n— >oo 

for every A; G N. 

As we already mentioned in the beginning of this section, the proof of existence of a 
minimizer of (3.2) uses an induction argument which depends on some mild qualitative 
property of the minimizer (see Theorem 3.8 and Theorem 3.9). Thus, we will prove first 
a stronger result which states that any eventual solution of (3.2) is a bounded set. We 
introduce the following notion which turns out to be a powerful tool in the analysis of 
optimal sets (see [7] [TO] and [TO] for similar techniques) 

Definition 3.4. Let Q, C R d be of finite perimeter and finite Lebesgue measure. We say that 
$7 is a local energy subsolution with respect to the perimeter or simply energy 
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subsolution if there are constants e = e(O) > and k = fc(ft) > such that for each 
measurable set [/ C ft with the property 

\\wn - wu\\ L 2 < e, (3.4) 

where wq and wjj are the solutions of ( |2.5[ ), we have 

E(n) + kP(Q)<E(U) + kP(U). (3.5) 

Proposition 3.5. Let f be a function satisfying the assumptions (fl), (/2) and (f3). Then 
any solution ft C M. d of (3.2) is an energy subsolution. 

Proof. Let U C ft and t = (P(n)/P(U)) 1/{d ~ l) . Suppose that t > 1, i.e. P(U) < P(ft). By 
the optimality of ft, properties (/2), (/3), the trivial scaling properties of the eigenvalues 
and of the perimeter and the monotonicty of eigenvalues with respect to set inclusion, we 
obtain 

o < / (x kl (tu) ,...,x kp (tuj) - f (\ kl (ft) , . . . , X kp (ft)) 
= / (\ kl (tU) ,...,X kp (tU)) - f (\ kl (17), ... , X kp (Uj) 

+ f (x kl (to, ... , x kp (t/)) - / (x kl (ft) , . . . , x kp (ft)) 
< ^(P(n))-^ fe^^j -P(n) 

+ L^(A fc ,(t/)-A fc ,(ft) 
i=i 

where L is the (local) Lipschitz constant of / and a is the constant from (/3). Using the 

2 

concavity of the function z i— > z d ~ x if c? > 3, or the fact that P(U) < P(ft) if d = 2, we can 
bound 

p(U)^i - P(ft)^r < C(ft) (P(t/) - P(ft)) . 



By Proposition 2.5 we have the estimate 

0<X kj (U)-X kj (Q)<C\\wn-wu\\l 2 , Vj = l,...,p. (3.6) 

Thus, there are constants A(ft) > and e(ft) > such that for each U C ft with the 
property 

P(t/)<P(ft), \\w a -wu\\^<e, 

we have 

v v 

^A fe .(ft)+AP(ft)<^A fc .(t/)+AP(t7). (3.7) 

3=1 3=1 

Moreover, thanks to the monotonicity of the eigenavlues by set inclusion, the above in- 
equality trivially holds also if P(ft) < P(U). 

We are now in the position to prove that any solution ft of (3.2 ) is an energy subsolution. 
Indeed, by Lemma 4.1], for every U C ft we have 



A fei (*7) - A*. (ft) < CX kj (U) (E{U) - E(Q)) 
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which together with (3.7) (which holds for every (/ C U satisfying (3.4)) and the mono- 
tonicity of the eigenavalues by set inclusion, give 

A (P(O) - P(U)) < A fc? (U) - X kj (Q) < CX k . (U) (E(U) - E(Q)) , 



where C is a constant depending on £1 and k p . Using again (3.6), we see that all the A&. (£/), 
for j = 1, . . . ,p, remain bounded by a constant depending only on f2 as E(U) — E(JY) < e, 
concluding the proof. □ 



Remark 3.6. In the case p = 1 and f(x) = x, the proof of Proposition 3.7 can be simplified 
due to the following fact: There exists a positive constant A > such that any solution 
ft C M. d of (3.1) is also a solution of the problem 



mm 



|A fc (0) + AP(O) : ncM d , Q measurable, |0| < +oo| . 



(3.8) 



In order to prove that, we first notice that, by the scaling properties of A& and of the 
perimeter, Q, C M rf is a solution of (3.8), if and only if, the following two conditions are 
satisfied: 

(1) n is a solution of ( pTT] ) with c = P(Q). 

(2) The function F(t) = Afc(ifi) + AP(tQ), defined on the positive real numbers, achieves 
its minimum in t = 1. 



Thus, if is a solution of (3.1 ), then it is sufficient to choose A > such that the derivative 

F'{t) = -2X k (n)t~ 3 + (d- l)AP{tt)t d - 2 , 

vanishes in t = 1, i.e. 

A - 2A ^) (39 ) 

Lemma 3.7. Let Q be an energy subsolution. Then f2 is a bounded set. 
Proof. For each t G M, we set 



H t = {x G R d : xi= t}, H+ = {x G R d : x\ > t}, H~ = {x G R d : x 1 < t} 



(3.10) 



We prove that there is some t G K such that n f2| = 0. For sake of simplicity, set 
w := and M = ||io||ioo. For any t G M. consider the function 



v t {xi, ...,x d ) 



M , xx < t - y/M, 

\{2M - (xi - t + V2M) 2 ^j ,t-V2M<x 1 <t, 
, t < x x . 



(3.11) 



Consider the set Hj = J!fl Hj~ , obtained "cutting" f2 with an hyperplane, and the 
function wt = w A vt G -f/g(fit). We recall that u>n t is the orthogonal projection of w on 
LiQ^t) with respect to the Hq(TI) scalar productj^j hence 

\\w nt -w\\ L 2 < C(n)\\Vw-Vw nt \\ L 2 < C(Q)||V(w - w t )\\ L2 , 



1 We recall that, thanks to the Poincare inequality, Hq(Q) is an Hilbert space with the scalar product 
given by 



9 



for some constant C(0) depending on f2. Thus, for t big enough, we have E(flt)—E(Q) < e. 
Hence, we can use fi 4 as a competitor in (13 . 51) . We thus get the inequality: 



j{w) + kP(n) = E(n) + kP(n) < E(n t ) + kP(n t ) < j( Wt ) + kP(n t ), 



where the functional J is defined in (2.6). Hence we get 

\ I \Vw\ 2 dx- I wdx + kP(Q)<- I \Vw t \ 2 dx- [ w t dx + kP(n t ). 
2 Jn Jn 2 J nt J nt 

Notice that wt = on and wt = w on H~ Setting t- = t — \/2M, and using the 



inequality 



we obtain 



\b\- 



< a ■ (a — b) V a, b G 



J/ \Vw\ 2 dx + k(P(n) - P(fl t )) < I [ \Vw t \ 2 -\Vw\ 2 dx+ [ (w-w r )dx 



< 



/ Vwt ■ V(wt — w) dx + I 

J{t-<xi<t} J 



(w — wt) dx 



{t-<X!<t} 

H t OXl 



'2M / wdU^ 1 + 



Vvt -V(w — vt)+dx + / (w — vt)+dx 



w dx 



w dx 



where for a generic u G if 1 (IR rf ), with u + := sup{w, 0} we indicate the positive part of u. 
Using again the boundedness of w, we get 



k(p(n, Hf) - p(h+, n)) < \p2M*i 2 u d - x {K t nn) + M|o dh+\. 

On the other hand, by the isoperimetric inequality, for almost every t we have 

\n n Hf\^ < c d p(n n h+) = c d (% d ~ 1 (H t n o) + p(q, 



Putting together (3.12) and (3.13) we obtain 



t i > 



(3.12) 



(3.13) 



(3-14) 



where Ci is some constant depending on the dimension d, the constant k and the norm M. 
Setting = \Qn Hf\, we have that c/>(i) as t -> +oo and = -H d ~ l (H t n O). 
Chosing T = T(ft) such that 

Ci0(i) < ^(i)T 1 Vt>T, 



equation (3.14) gives 



< -2C 1 0(t) 1 - 1 / d Vt > T, 



which implies that 4>(t) vanishes for some t G M. Repeating this argument in any direction, 
we obtain that 0, is bounded. □ 
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We are now in position to prove the existence of an optimal set for (3.2). We first prove 
the result for the problem (3.1). This is just a particular case of (3.2), but it will be the 
first step of the proof of the more general Theorem 3.9, which is based on the same idea. 

Theorem 3.8. For any k 6 N, there exists a solution f2 C M. d of (3.1). Moreover, any 
solution of (3.1) is a bounded set. 

Proof. Without loss of generality we assume that c = 1 . We prove the theorem by induction 
on k G N. For k = 1 the existence holds, since by a standard symmetrization argument, we 
have that the optimal set is a ball of perimeter 1 . 

Let k > 1 and let £l n be a minimizing sequence for (3.1). We note that clearly the 
perimeters are bounded and the sets have finite measures, hence the assumptions of Theo- 
rem 



3.1 are satisfied. Moreover by the monotonicity of the eigenvalues 
limsupAi(0„) < limsup Ajt(f2 n ) < +oo, 



hence by Remark 3.2 we have only two possibilities: 

(i) Compactness: Since y n + Q n is also a minimizing sequence, we have that the limit Q 
is such that 

P(Q) < 1 and X k (Q) < hminf A fc (fi n ). 

n— i>oo 

Thus Q is a solution of (|3.1 



(ii) Dichotomy: By the dichotomy case of Theorem 3.1 and a scaling argument, we can 
construct a sequence of sets 

A n U B n , 

satisfying P(Q n ) = 1, d(A n ,B n ) — > oo and which is still a minimizing sequence for 
(3.1). Without loss of generality and up to extracting a subsequenca^J if necessary, 
we may assume that X k (A n U B n ) = Xi(A n ) for some I < k. We note that since Q, n 
is minimizing, we must have I < k. Indeed, if this is not the case, i.e. I = k, the 



sequence A n is such that X k (A n ) = Aj.(O n ) but on the other hand, by Theorem 

limsupP(A n ) < 1 - liminf P(B n ) < 1, 



3.1 



where the strict inequality is due to point (c) of the dichotomy case in Theorem 3.1 
and the isoperimetric inequality 

liminf P(B n ) > C d lim inf | B n \ ^ > 0. 

n— >oo n— >co 

Thus an appropriate rescaling of A n would lead to a strictly better minimizing se- 
quence, which is impossible. Hence I < k, taking the biggest possible I, we can assume 
A;(^4 n ) < A; + i(^4 n ). Since the spectrum of the Dirichlet Laplacian of Q n is given by 
the union of the spectrum of the the Laplacian on A n and B n , this assumption implies 



max • 



.{Xi(A n ),X k ^i(B n )} = Xi(A n ) = A fc (fi„). 

Up to a extract a subsequence, we suppose that the following limits exist: 
lim P(A n ), lim P(B n ), lim Xi(A n ), 



(3.15) 



lim X k -i(B r , 

n— >oo 



2 We recall that if Q. = A U B with dist(,4, B) > 0, Hq(Q) = Hq(A) Hq(B) hence the spectrum of the 
Dirichlet Laplacian of 57 is given by the union of the spectrum of the Dirichlet Laplacian of A and of B. 
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Let O* and be the optimal sets for ( |3.1| ) for A; with constraint ca = hmn._j.oo P(A n ) 
and Afc_; with constraint cb = limn^oo P{B n ), respectively. Since, thanks to Lemma 



3.7, they are bounded, we may suppose that up to translation, they are at positive 



lim P{A n ) + lim P(B n ) < 1 

n— ¥00 n—^oo 



distance. Then 

and, by the choice of fij and 

A^) < lim Xi(A n ), Xk-ii^U) < lim X k -i(B.. 



As a consequence, thanks to (3.15), 



A fc (nj" U = max (Aj(n*), Ajfe_,(J_3-_ £ )} 

< lim max{A;(A n ), A fc _/(-B„)} 

= liminf Afc(fJ n ). 



Hence f2* U fij^ is a solution of (3.1). 



□ 



We now prove the existence in the general case of a spectral functional of the form 

F(n) = f(\ 1 (n),...,\ k m 



Theorem 3.9. Suppose that f satisfies (/l), (/2) and (/3). T/ien there exists a solution 
n C R d of dO). Moreover, every solution of (|3.2h is a bounded set. 



Proof. As in the proof of Theorem 3.8 we proceed by induction, this time on the number 
of variables p. If p = 1, then thanks to the monotonicity of /, any solution of (3.1) is also 



a solution of (3.2) and so we have the claim by Theorem |3.8| 

Consider now the functional F(£l) = /(Afe 1 (0), . . . , Afcp (O)) and let O n be a minimizing 
sequence. By Theorem 3.1 and using (/l) and Remark 3.2 there are two possible behaviours 
for the sequence Q. n : compactness and dichotomy. 

If the compactness occurs, we immediately obtain the existence of an optimal set. Oth- 
erwise, in the dichotomy case, we may suppose that Q n = A n U B n , where the Lebesgue 
measure of A n and B n is uniformly bounded from below and dist(^4 n , B n ) —> oo. Moreover, 
up to a extracting a subsequence, we may suppose that there is some 1 < I < p and two 
sets of natural numbers 



1 < a\ < ■ ■ ■ < ai, 
such that for every n € N, we have 

jA ai (Ai)j • • • > X ai (A n ), X/3 l+1 (B n ), . 



i < A+i <■■■<#,, 



^/3 p (Bn)\ = jAfc^^n), . . . , Afc p (^n)} • 



Indeed, if all the eigenvalues of O n are realized by, say, A n arguing as in the proof of 
Theorem |3.8| we can construct a strictly better minimizing sequence. Moreover, without 
loss of generality we may suppose that 

X ai {A n ) = A fc .(O n ), Vi = 1, ... ,1; Xp.{B n ) = X kj (Q n ),Vj = 1 + 1,..., p. 

We can also assume that the sequences X ai (A n ) and Xp.(B n ) converge as n — > oo. By a 
scaling argument as in Theorem 3.8 we also have that without loss of generality P(A n ) = c a 
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and P(B n ) = cp, where c a and cp are fixed positive constants. Let f a : M. 1 
restriction of / to the Z-dimensional hyperplane 



be the 



\x G R p : Xj = lim Xp. {B n ), j = I + 1, . . . ,p\ . 
Since I < p, by the inductive assumption, there is a solution A* of the problem 

min|/ a (l\ ai (A), . . . ,X ai (A)) : A C M d , A measurable, P(A) = c a , \A\ < ooj . (3.16) 
Since / is locally Lipschitz, we have 

liminfn^oo / (X ai (A n ), . . . , X ai (A n ), Xp l+1 (B n ), . . . , Xp p (B n ) 



liminf^oo / f X ai (A n ), . . . , X Ql (A n ), limm^oo A^ J+1 (B m ), . . . jUm^oo Xp p (B n 



and thus the minimum in (3.16) is smaller than the infimum in (3.2). Moreover, A* is 



bounded and so, we may suppose that dist(^4*, B n ) > 0, for all n S N. Thus, again by the 



Lipschitz condition on /, the sequence A* U B n is minimizing for (3.2). 



Let now fp :M. P 1 — > R be the restriction of / to the (p — /)-dimensional hyperplane 
{xef: Xi = X ai (A*), i = l,...,l}. 
Let B* be a solution of the problem 



mm 



{fp(Xp l+1 (B),...,Xp p (B)\ : BcR d , B measurable, P{B) = cp, \B\ < oo} 



(3.17) 

We have that the minimum in (3.17) is smaller than the minimum in (3.16) and so than 



that in (3.2). On the other hand, since both A* and B* are bounded and the functionals 



we consider are translation invariant, we may suppose that dist(A* , B*) > 0. Thus the set 
Q* := A* U B* is a solution of ((331). □ 



4 Existence of an open solution 



In this section we study the shape optimization problems (1.7) and (1.8). In Theorem 



4.8 we prove that there exist solutions of (1.7) and (1.8) and we study their relation with 



the corresponding solutions of (3.1) and (3.2). Before we prove the theorem we need some 



preliminary results concerning the sets which, in some generalized sense, have positive mean 
curvature. 

Definition 4.1. We say that the measurable set 0, is a perimeter supersolution if it 

has finite Lebesgue measure, finite perimeter and satisfies the following condition: 



P(D) < P(n), for each SlDll. 



(4.1) 



Remark 4.2. Let f2 be an open set with boundary dfl of class C 2 . If 17 is a perimeter 
supersolution, then it has positive mean curvature with respect to the exterior normal 



vector field on dft. Lemma 5.3 below shows that, even if it is less regular, it has positive 



mean curvature in the viscosity sense. 

The following simple Lemma will play a crucial role in the sequel: 
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Lemma 4.3. Suppose that f : W — > R satisfies conditions (/l), (/2) and (/3) and i/iai 
£7 C M d is a solution of (3.2). T/ien n is a perimeter supersolution. 



Proof. Suppose, by contradiction, that D is such that P(O) < P(O) and set 

/ ~ x l/(d-l) 

t = (p(n)/p(n)j >i. 

Then, for any k G N, we have 

A fe (iO) < A fc (0) < A fc (J2). 

On the other hand P(tQ) = P(O) and so, by the optimality of Q and the strict monotonicity 
of /, (/3), we have 

o </(A fcl (tn),...,Afc p (tn)) -/(AfcCn),...,^^)) 
</(A fcl (n),...,A fci) (n))-/(A fcx (n),...,A* p (n)) <o, 

which is a contradiction. □ 

The following result is classical (see, for instance, [23], |30|, Theorem 16.14]) and so we 
only sketch the proof. 

Lemma 4.4. Let £7 C M. d be a perimeter supersolution. Then there exists a positive constant 
c, depending only on the dimension d, such that for every x G M. d , one of the following 
situations occur: 

(a) there is some ball B r (x) with r > such that B r (x) C SI a.e., 

(b) for each ball B r (x) C R d , we have \B r {x) n 9, c \ > c\B r \. 

Proof. Let x G M. d . Suppose that there is no r > such that B r (x) C 0,. We will prove 



that (6) holds. Using the condition (4.1 ) for O = fiU B r (x) we get that for almost every r, 

p(n, B r {x)) < n d - l {dB r {x) n n c ). 

Applying the isoperimetric inequality to B r (x) \ f2, we obtain 

\B r (x) \ ni 1 - 1 ^ < C d (P(SI, B r (x)) + U d - l (dB r {x) n o c )) 

(4.2) 

< 2C d n d - 1 (dB r (x) nn c ). 

Consider the function <j>{r) = \B r (x)\Sl\. Note that (f)(0) = and cp'(r) = U d ~ l {dB r (x)C\Q) 
and so, by (4.2 ), 

c< 4-(4>(r)W 



dr 

which after integration gives (b). □ 

Definition 4.5. If ft C M. d is a set if finite Lebesgue measure and if there is a constant 
c > such that for each point x G M. d one of the conditions (a) and (b), from Lemma 
holds, then we say that SI satisfies an exterior density estimate. 



4.4 
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Proposition 4.6. Let C M. d be a set of finite Lebesgue measure satisfying an exterior 
density estimate. Then there are positive constants C and ft such that, for each x £ M d 
with the property that \B r {x) n £l c \ > 0, for every r > 0, we have 

\\wn\\L™(B r {x)) < rf3 \\wn\\L°°(R d )i f° r each r >0. (4.3) 

In particular, if Q is a perimeter supersolution, then the above conclusion holds. 

Proof. Let x 6 M. d be such that that \B r (x) n Q c \ > 0, for every r > 0. Without loss of 
generality we can suppose that x = 0. Setting w := ioq, we have that Aw + 1 > in 
distributional sense on M. d . Thus, on each ball B r (y) the function 

u(x) = w(x) - ^(r 2 -\x - y\ 2 ), 

is subharmonic. By the mean value property 

r 2 1 

2d u d r d J Br{y) 



w{y) < — H -g / w(x)dx. (4.4) 



Let us define r n = 4 n . For any y S -B rn+1 (0), equation (4.4) implies 



r 2 If 

4d |B 2 r„ +1 (2/)| JB 2rn+1 (y) 



rl , inns^Cy)! 



< + 



4d |S 2r „ +1 (y)| 



w||l~(b 2 (B) ) 



(4.5) 



. r 2 / |n c na rw+1 (o)| ^ 

4" 2n / _ rf _\ „ , 

<^" + (l"2 C ) IMU«(B^(0)), 

where in the third inequality we have used the inclusion B rn+1 (0) C B2 Tn+1 {y) for every 
y E B rn+1 (0). Hence setting 

a n = IMU°°(Br n (0)), 

we have 

g— n 

a n +l < -jv + (1 - 2~ d c)a n , 
which easily implies a n < Cao4~ nl3 for some constants [3 and C depending only on c. This 



gives (4.3). □ 



Proposition 4.7. Let C M. d be a set of finite Lebesgue measure satisfying an external 
density estimate. Then the set 

!!i:= leT: 3 lim l — = 1 

[ r '^° \B r (x)\ 

is open and Hq(Q) = Hq(£Ii). In particular, if f2 is a perimeter supersolution, then J)i is 
open and Hq(Q) = ifg(Jli). 



15 



Proof. Thanks to Lemma 4.4 J7i is an open set. It remains to prove the equality between 
the Sobolev spaces. We first note that we have the equality 

H 1 (n) = H 1 ({w n >0}), 
where wq is as in Section [2J Indeed, by the definitions of Hq and Hq, we have 

H^{w a > 0}) c Hl{{w a > 0}) c H^(n), 

and so it is sufficient to prove that the inclusion Hq(Q) C Hq({wq > 0}) holds. For, it is 
enough to check that for every positive and bounded u G Hq(Q) we have u G Hq({wq > 0}). 
Reasoning as in [20, Proposition 3.1], for every u G Hq(Q) such that < u < 1 and every 
n G N, we consider the weak solution u n G Hq(Q) of the equation 

— Au n + nu n = nu. 

By the weak maximum principle, we have that u n < nwn a.e. and so, by |26[ Lemma 
3.3.30], u n < hwq quasi-everywher^J which implies that u n G Hq({wq > 0}). On the 
other hand, using u n — u G Hq(Q) as a test function, we have 



/ Vu n • V(w n — u) dx + n I \u n 
Jn Jn 



u\ 2 dx 



and so 



/ \V(u n — u)\ 2 dx + n / \u n — u\ 2 dx = — Vu • V(u n — u) dx < \\Vu\\ L 2\\V(u n — u)\\ L 2. 
Jn Jn Jn 

(4-6) 

By (4.6), we have that ||V(n n — u)\\ L 2 < ||Vit||^2 and \\u n — u\\ < n 1 ||Vn|| i 2. Thus u n 
converges to u weakly in Hq(Q) and, by the first equality in (4.6), the convergence is also 
strong in As a consequence u n converges pointwise to u outside a set of zero 

capacity (see [Ml Proposition 3.3.33]) and so we have u G ^({wj) > 0}). 

We now prove that Q± = {wn > 0} up to a set of zero capacity. Consider a ball B C Oi. 
By the weak maximum principle, uig < wn and so 

VLi C {wn > 0}. 

In order to prove the other inclusion, we first note that since wn G iif 1 (M rf ), then there is a 
set iV C W 1 of zero capacity such that for any xq G M d \ N there exists the limit 



lim — — — -r / wn dx =: wn(x ), 
r->0 \B r (x )\ J Br (x ) 

and wn coincides with wq again up to a set of zero capacity (see [2H Section 4.8]). By 



Proposition 4.6, wn = on M \ f2i which gives the other inclusion. □ 



Theorem 4.8. Let k\, . . . , k p G N and suppose that f : W — > M satisfies the assumptions 
(fl), (/2), and (/3). Then there exists a solution Q, C M d of (1.8). Moreover, we have 
that 



3 A property V is said to hold quasi-everywhere if 



cap^T 3 is false}) = 
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(i) every solution of (1.8) is a solution of (3.2); 



(ii) every solution of (3.2) is equivalent to a solution of (1.8). 



4.7 



Proof. Let f2 be a solution of (3.2) and let fii be as in Proposition 
open set with perimeter P{U) = 1. Using the inclusion Hq(U) C Hq(U), the optimality of 



Let U C 



be an 



Q and Proposition 4.7 we have that for any k G N, 



A fc (l7) > A fc (l7) > A fc (n) = A fc (ni), 



(4.7) 



i.e. fii is a solution of (1.8), which proves the existence part and claim (ii). Suppose now 



that U is an optimal set for (1.8). Then all the equalities in (4.7) must be equalities and 
so, we have claim (i). □ 



5 Regularity of the free boundary 

In this section we study the regularity of the reduced boundary <9f2 of any solution Q of 
(1.8). The goal is to show that solutions of (1.8) are quasi- minimizers for the perimeter 
and then to apply some classical regularity results. In order to prove the quasi- minimality 

P(n,B r ) < P(n,B r ) + Cr d , VOAficB r (4 



of the optimal set $7, we first note that, thanks to Propositions |4.6| and |4,7[ wq is continuous. 
Moreover since f2 is a perimeter supersolution we can show that dist(x, O c ) is super harmonic 
in £1 in the viscosity sense. Hence we can apply the maximum principle in order to prove 
that wn is Lipschitz. An estimate on the perimeter of the variations of an energy subsolution 
will finally give the desired quasi-minimality property of f2. The following lemma is classical 
and so, we only give a reference for the proof. 

Lemma 5.1. Let u G H 1 (B r ) be such that —Au = f G L°°(B r ) in a weak sense in the ball 
B r . Then we have 

l|Vn|| L oo (jBr/2 ) < C d ||/|| L oo (jBr) + —\\u\\ L oa( Br y (5.1) 

Proof. See [22]. □ 

Proposition 5.2. Let C M. d be a perimeter supersolution. Then wn : R d ^ R is Holder 
continuous and 

\w Q (x) - w Q (y)\ < C\x - yf , (5.2) 
where (3 is the constant from Proposition |4.6[ 



Proof. Thanks to Proposition (4.7), up to a set of capacity zero, we can assume that 
J7i is open and that wq is the classical solution, with Dirichlet boundary conditions, of 
—Awq = 1 in Sl]_. Consider two distinct points x,y G R d . In case both x and y belong to 
f2f, the estimate (5.2) is trivial. Let us assume that x G and let xo G dQ\ be such that 

\x — xq\ = dist(x, dil{). 

We distinguish two cases: 
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Suppose that y G M. d is such that 

2|x — y\ > dist(x, 90i). 



Hence x,y G -B^-j,] (xo) and by Proposition 4.6, we have that 

< C|x - and wn(y) < C\x - yf. 

Thus we obtain 

\wn{x)-w n (y)\<2C\x-yf. (5.3) 
Assume that y G M d is such that 

2|x — y| < dist(x, 



Applying Lemma 5.1 to wq in -Bdist^dQi)^) ^ we obtain 



l|Vton|U» (BdiBt(B|ani)/a(a!)) < dist(a . >ani ) ^ Cddistfodfii)" , (5.4) 

which, since /3 < 1, together with our assumption and the mean value formula implies 
\wq(x) - w n (y)\ < C d dist(x, dQi)P~ l \x -y\<\x- yf. 

□ 

In the following Lemma we show that a perimeter supersolution has positive mean 
curvature in the viscosity sense. This is done showing that the function d(x, O c ) is super 
harmonic in SI in the viscosity sense (see [T7] for a nice account of theory of viscosity 
solutions). In case d£l is smooth this easily implies that the mean curvature of d£l, computed 
with respect to the exterior normal, is positive (see for instance [22, Section 14.6]). A similar 
observation already appeared in |18j . in the study of the regularity of minimal surfaces, and 
in [27, 32 1, in the study of free boundary type problems. 

We say that (p touches da from below at xq if 

dn(xo) - <p(xo) = min {cfo - ip\. 

n 

Lemma 5.3. Let £1 C M. d be a perimeter supersolution. Consider the function cfo(x) = 
d(x, O c ). Then for each ip G C£°(f2) ; touching dn from below at x$ G ft, we have A<^(xo) < 
0. 

Proof. Suppose, by contradiction, that there are point xq G £1 and a function (p G C^°(J1) 
touching d^i from below at xq for which A(p(xq) > 0. Up to a vertical translation, we can 
assume that ip(xo) = dfi(xo) > 0. Furthermore, by considering a a regularized version of 
the function ip(x) = (<p(x) — \x — xo| 4 ) + , we can also suppose that (p(x) < dn(x), for every 
x G £l different from xq. 

Let yo G dQ be such that ofo(xo) = \xo — yo\, then 

V^o) = (5.5) 
\xo - Vol 

In order to prove this last equality, we first notice that, since (j) is smooth, the inequality 

(p(x) - cp(x ) < dn(x) - d n (x ) < \x - x \, (5.6) 
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gives |V</?(xo)| < 1. Moreover, defining xt ■= txo + (1 — i)yo> we have 

<ki(xt) = \xt - yo| = tdn(x ), 

hence 

~ , , x -y ip{x t ) - (p{x Q ) dn(x t ) - dn{x ) 

— Viplxo) ■ | r = iim — : : — < iim : : 

Fo-yo F* -a? o t -^ 1 Fi-^o 



which, together with (5.6), proves (5.5). 

Let us now set h := <p(xo) = cfo(xo) an d choose a system of coordinates such that 
xq = and the unit vector e d is parallel to xq — uq. Since 7^ 0, by the Implicit Function 
Theorem, there is a (d — l)-dimensional ball Bf~ l C M rf and a function G C°°{Bf~ 1 ) such 
that {99 = h} is the graph of (f) over B^ 1 , i.e. 

= /i} n (s*" 1 x (-r, r)) = {x d = 0(xi, . . . , Xd-j)}. (5.7) 

Since dn > tp with equality only at xq = 0, we have 

W > h} C {d n > h} and = h} n {(fo = /i} = {0}, (5.8) 

which implies that is a (strict) local minimum of the function 

(xi, . . . , x d -x) (->• x 2 H h x 2 ,^ + (0- /i) 2 . 

Hence J^-(0) = • • • = g§~(0) = 0. On the other hand, since 

V?(xi, . . . , x d _i, 0(xi, . . . , acd_i)) = 0, 

we get, denoting with the subscripts the partial derivatives, 

Lfj + <j)j(p d = 0, 

fjj + tyjVjd + 4>jjfd + 4> 2 jUdd = 0, 

for each j = 1, . . . , d— 1, and thus we obtain (fjj(0) + <f>jj(0)<p d (0) = 0- By the contradiction 
assumption 

d d-l d-l 

< <PiM = Vdd(0) - ^(o) 52 0ii(O) < -Mo) 52 

where the last inequality is due to 

(Q) = Hm <p(te d ) + f(-te d ) - 2y?(0) < ^ d n (te d ) + d n (-te d ) - 2d n (0) < Q 
t 2 ~ i-+0 i 2 



Since, by (5.5), we have ip d (0) = 1, we deduce that A<^>(0) < 0. 

Let ds : T + — >• M be the distance to the surface S = {x d = <j)(xx, . . . , x<j_i)}. i.e. 
ds(x) = d(x,S), where T is a tubular neighbourhood of S 1 and T + = {x d > <fi}. Then 
d s £C°°{T + US), 

Arguing as above, we see that Ads(0) = —A<f>(0) > and so, Ads > 0, in a neighbourhood 
of in T+ U S. 
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S c = {ip = h] - (h + e)e d 



d(l = {d n = 0} 



Figure 1: Proof of Lemma |5.3| applying the Divergence Theorem to the grey region, we 
obtain a contradiction to the minimality of f2 if Ay? > 0. 



By (5.8) we see that for r small enough, there is some e > such that 

{h < d Q < h + e} n {<p > h} C B r . 

If we define the set 

fie :=fiU ({ip >h}-(h + e)e d ) , 
then fi e \ Q C B r (—(h-\- e)e<j). Denoting with d 6 the distance from 

S € = {tp = h}-(h + e)e d , 

we see that Ad e > in B r (—(h + e)e^), since d t (x) = ds(x + (h + e)ed). Hence, by the 
Divergence Theorem, and recalling that on S t , Vc? e = — where is the exterior 
normal to O e , we have 

< / Ad e dx = - [ Vd e - v a dU d - x - [ dU d ~ x 



n £ \n Jn £ nen Jn c ndn e (5.9) 

< u d " x (n e n an) - n^ 1 (n c n dn c ) , 

contradicting the perimeter minimality of Q with respect to outer variations (see Figure 
0. □ 

We are now in position to prove the Lipschitz continuity of wn using do, as a barrier 
(see |22l Chapter 14] for similar proofs in the smooth case). 

Proposition 5.4. Suppose that the open set Q C M. d is a perimeter supersolution. Then 
the energy function wq : M rf — > M., defined as zero on £l c , is Lipschitz continuous. 

Proof. For sake of simplicity, we set w = wn and || • ||oo = || • ||L°°(n)- Let c > 2||w||oo and 
consider the function 

h(t) = ct-t 2 . (5.10) 

We claim that 

w(x) < h(d n (x)) VxeU. (5.11) 
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Suppose this is not the case. Since both functions vanish on dQ, there exists xq G f2 such 
that 

h(dn(x )) - w(x ) = min {h(d n ) - w\, 

n 

that is the function 93 := h~ 1 (w) touches efo from below. By our choice of c the function 
h is invertible on the range of w. Moreover, since wq(xq) > 0, the inverse function is also 



smooth in a neighborhood of xq. By Lemma 5.3 



A^(x ) < 0. 



Hence, the chain rule and the definition of h, (5.10), imply 

Aw(x ) =A(ho<p) (x ) = h"(<p(x ))\X7<p(x )\ 2 + ti{tp(x ))A<p(x ) < -2|V^(a:o)| 2 = -2, 

where we have also taken into account that, since (p touches dn from below at xq, equation 
5.5) implies the |V^(xq)| = 1. Since — Aw = 1 the above equation cannot hold, hence 



5.11) holds true. Now equations (5.11) and (5.10), imply 



w(x) < h(dfi(x)) < cdn(x) \/ x £ S7. 



Arguing as in the proof of Proposition 5.2 we conclude that w is Lipschitz. 



□ 



In order to prove that any solution of (1.7) is a quasi-minimizer of the perimeter we 
need to use both the facts that it is a perimeter supersolution (see Lemma 4.3) and energy 
subsolution (see Proposition 3.5). In the next lemma, we will obtain some local information 
on dQ using an inner variation of f2. This lemma in combination with Lipschitz continuity 
of f2 will give the regularity of the boundary of (see Theorem 5.6 below). The argument 
below is classical and similar results can be found in [2], [7] and [10J. 

Lemma 5.5. Let Q, C R rf be an energy subsolution and let w = wq. Then for each 
B t (xq) C M rf and each Q, C f2 such that 17 \ C B r {xo), we have the following inequality: 



\Vw\ 2 dx + k(P(n) - P(O)) 



Br(x ) 



< 



wdx + Cd[r + 



B r (x ) 



M\L°°{B 2 r{x )) 



u 



wdH d -\ (5.12) 



8B r (x ) 



Proof. Without loss of generality, we can suppose that xq = 0. We will denote with B r the 
ball of radius r centered in and with A r the annulus B<i r \ B r . 
Let ij) : A\ — > M + be the solution of the equation: 

Aip = on A±, V = on dB±, ip = 1 on 8B2- 

We can also give the explicit form of tp, but for our purposes, it is enough to know that tjj 
is bounded and positive. 

With <p : A\ — > M + we denote the solution of the equation: 



: 1 on 4i, = on dB±, <ft = on dB2- 
For an arbitrary r > 0, a > and k > 0, we have that the solution of the equation 

— Av = 1 on A r , v = on dB±, v = a on 8B2, (5.13) 
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is given by 

and its gradient by 



v(x) = r 2 (p(x/r) + aip(x/r), 



OL 

Vv(x) = rV(j)(x/r) -\ — Vip(x/r). 

r 



(5.14) 



(5.15) 



Consider the solution v of the equation 5.13 with a = ||w||^oo(5 2r ). We work with the 

perturbation w r = w\bs + w A fl_B 2r . Observe that, by the choice of a, w r G Hq(Q,). 
Moreover, a simple computation gives 

\E(tt) - E(Q)\ ->■ uniformly as r ->• 0. 

Being an energy subsolution, we get, for r small enough, 



Vw\ z dx- / u>(x) dx + feP(O) < 



1 



X7w r \ 2 dx — / w r (x) c?x + kP(£l). 



Since uv = in B r and wv = w in (-E>2r) c 5 we obtain with the same computations of Lemma 
X71 



|Vw| 2 + k(p(n) - p(n)) < 



i 



I VWr I 



I Viol dx + 



B 2 r 



(w — uv) <ix 



< 



V« • V((u> — w) V 0) dx + / (w — w r ) dx 

r J B 2 r 

dB r on 



((w — v) V 0) dx + / (w — w r ) dx 



dB r dn 



w dx 



13 r 



at , 



< (r||V0Hoo + -[|V^[| c 



w dH 



d-l 



+ 



dB r 



B, 



B 2r 

w dx, 
(5.16) 



where the last inequality is due to (5.15). Recalling that a 
<p2b. 



|w||£oo(£ 2r ), we obtain 
□ 



Theorem 5.6. Let Q C M d be a set of finite Lebesgue measure and finite perimeter. If Q 
is an energy subsolution and a perimeter supers olution, then Q is a bounded open set and 
its boundary is C 1,a for every a £ (0, 1) outside a closed set of dimension d — 8. 

Proof. First notice that, by Lemma |3.7[ f2 is bounded. Moreover, since f2 is a perimeter 
supersolution, we can apply Proposition 4.7 and Proposition 5.4, obtaining that f2 is an 
open set and the energy function w := wn is Lipschitz. 
We now divide the proof in two steps. 



Step 1. Let xq £ dCl and let B r (xo) be a ball of radius less than 1. By Lemma 5.5 for 



each f2 C 0, such that £1AQ C B r (xo), equation (5.12) implies (recall that r < 1) 

wdn^ 1 



k(p(n) - p(n)) < 



w dx + Ca r + 



B r (xo) 



\w\\L°°(B 2 r(xo)) 



dB r (xo) 



(5.17) 



< C d \\w\\ L oo( B2r ( XQ ^r X > 

where Cd is a dimensional constant. Now since w is Lipschitz and vanishes on <9f2, we have 
ll' u; llL 00 (B2r(^o)) — hence equation (5.17), implies 

P(0, fl r (x )) < P(fi, B r (x )) + Cr d , (5.18) 
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where C depends on the dimension, the constant k on (3.5) and the Lipschitz constant of 
w (which, in turn, depends only on the data of the problem). Moreover, by the perimeter 
subminimality, equation (5.18) clearly holds true also for outer variations. Splitting every 
variation in an outer and inner variations (see for instance |19t Proposition 1.2]), we obtain 

P(Sl,B r ) < P(n,B r ) + Cr d , VflAficB r (4 

Hence is a almost-minimizer for the perimeter in the sense of |35U36| (see also pQ). From 
this it follows that dSl is a C 1,a manifold, outside a closed singular set £ of dimension 
(d-8), for every a E (0,1/2). 

• Step 2. We want to improve the exponent of Holder continuity of the normal of dSl in the 
regular (i.e. non-singular) points of the boundary. For this notice that, for every regular 
point xo £ dSl, there exists a radius r such that dSl can be represented by the graph of a 
C 1 function <j> in B r (xo), that is, up to a rotation of coordinates 

SI n B r (x ) = {x d > (p(xi, Xd-i)} n B r (x ). 

For every T £ Cl(B r (xo)',^- ) such that T • vq < and t is sufficiently small, we consider 
the local variation 

Sit = (Id + tT)(Sl) C SI. 
By the energy subminimality we obtain 



k(P(Sl) - P(Sl t )) < E(Sl t ) - E(Sl). 



(5.19) 



Since T is supported in B r and dSl n B r is C , we can perform the same computations as 
in [56l Chapter 5], to obtain that 



E(Sl t ) -E(Sl) =-t f 
Jd 

Moreover, see for instance \30\ Theorem 17.5 



/ 




1 dnnBr 


dv 



T ■ i/fi dU d - 1 + o(t) 



P(Sl t ) = P(Sl) + t ( div m T dU d ~ x + o{t) 
J dnnBr 



(5.20) 



(5.21) 



where divgnT is the tangential divergence of T. Plugging (|5.20|) and (|5.21|) in (|5.19|), a 

S3 (in 

1 dwn 



standard computation (see |3Q|. Theorem 11.8]), gives (in the distributional sense 

2 



div 



< 



dv 



where the last inequality is due to the Lipschitz continuity of wq. Moreover applying (5.21 ) 
to outer variations of SI (i.e. to variations such that T • vq > 0) we get 



div 



> 0. 



In conclusion is a C 1 function satisfying 

div ( - ^ ) £ l- 

and classical elliptic regularity gives (p G C 1,a , for every a £ (0, 1). 



□ 
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Remark 5.7. In some particular cases the conclusion of the above Theorem can be strength- 



ened. For instance, if one knows that, for a solution of (1.7), the kth eigenvalue is simple, 
then one can write the Euler Lagrange equation for the minimum domain which reads, in 
a neighbourhood of a regular point, as 



div 





dw k 


)■ 


dv 



where Wk is the kth. eigenfunction. In this case a standard bootstrap argument gives that 
the regular part of d£l is an analytic manifold. 









"e~ 


fi 











2g 



Figure 2: The variation in the proof of Proposition 5.8 



Proposition 5.8. Let f : W -) 

ficM^ be a bounded open set, solution of (1.8). Then Q is connected. 



be a function satisfying (fl), (/2) and (/3) and let 



Proof. We first prove the result in dimension d < 7, in which case the singular set of the 



boundary dfl is empty. We first note that, since Q is a solution of (1.8), it has a finite 
number (at most max{ki, . . . , k p }) of connected components. Suppose, by contradiction, 
that there are at least two connected components of Q. If we take one of them and translate 



it until it touches one of the others, then we obtain a set Q which is still a solution of (1.8). 
Using the regularity of the contact point for the two connected components, it is easy to 
construct an outer variation of which decreases the perimeter (see Figure 0). In fact, 
assuming that the contact point is the origin, up to a rotation of the coordinate axes, we 
can find a small cylinder C r and two C 1,a functions g\ and g 2 such that 



and 



51 (0) = g 2 (0) = |V 5 i(0)| = |V 52 (0)| = 0, 



Q, c nC r = {gx{xx,...,x d -\) <x d < g 2 (xi,...,x d - 1 )} nc r 



(5.22) 



Now, for q < r, consider the set Q e := Q, U C g D It is easy see that, thanks to (5.22) 
and the C 1,a regularity of g\ and g%, 

p(n g ) - p(o) < c a0 d - l+a - c d6 d - 1 < o, 

for g small enough, which contradicts the minimality of Q 



4 Another way to conclude is to notice that for fl the origin is not a regular point, a contradiction with 



Theorem 5.6 
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We now consider the case d > 8. In this case the singular set may be non-empty and 
so, in order to perform the operation described above, we need to be sure that the contact 
point is not singular. 

Suppose, by contradiction, that the optimal set O is disconnected, i.e. there exist two 
non-empty open sets i,5cSJ such that A U B = and A n B = 0. We have 

Mi 



dAudB c dn = d M n, 

where the last inequality follows by classical density estimates. By Federer's criterion [30|. 
Theorem 16.2], A and B have finite perimeter. Arguing as in [3, Theorem 2, Section 4], we 
deduce that P(O) = P{A) + P(B). 

Since both A and B are bounded, there is some xq G M d such that dist(A, xq + B) > 0. 



Then the set f2' = A U {xq + B) is also a solution of (1.8). Let x G cL4 and y G <9(xo + -B) 
be such that |x — y\ = dist(^4,xo + B). Since the ball with center (x + y)/2 and radius 
\x — y\/2 does not intersect Q 1 , we have that in both x and y, £1' satisfies the exterior ball 
condition. Hence both x and y are regular point^J 



Consider now the set f2" = (— x + A) U (— y + xq + B). It is a solution of (1.8) and has at 
least two connected components, which meet in a point which is regular for both of them. 
Reasoning as in the case d < 7, we obtain a contradiction. □ 



6 Appendix: Proof of Theorem 3.1 



We apply the concentration compactness principle from [29] to the sequence of characteris- 
tics functions ln n - First notice that, being all the sets of finite measure, the isoperimetric 
inequality and the uniform bound on the perimeters ensure that 

sup \ Q n \ < C. 

nt=N 

Let us assume that 

limsup |f2 n | > 0. (6.1) 
In this case, up to subsequences, we can assume that 

/ := lim |$7 n | 

exists and that / G (0,+oo). Thanks to this can rescale all our sets in such a way that 
|ri n | = 1 still maintaining a uniform bound on the perimeters. 

As in [29] we have that, up to subsequences, the family of concentration functions 
Q n : R + M+, defined by 

Qn(r) = sup \B r (x) n fl n \, 



is pointwise converging to some monotone increasing function Q : M + — > M + . We now 
consider three different cases: 

(i) linv^oo Q(r) = 1. In this case, up to substitute Q n with Q n + x n for suitable x n G M. d , 
we have that for every e > there is some R > such that 

sup |O n \ B R \ < e. 



5 We denote with d M Q the essential boundary of ft, i.e. the complement to the set of density 1 points of 
Q and of fl c . 

6 This can be easily seen, since any tangent cone at these points is contained in an half-space and hence 
it has to coincide with it, see [331 Theorem 36.5] 
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Since the functions Wfi n are uniformly bounded in L°°(]R a! ) we infer that for every 
e > there is some R > such that 

sup / wn n dx < e. 

neN J B C R 

By the compact inclusion BV(R d ) ^ Ll oc (R d ) and fl' 1 (R d ) ^ Lf oc (]R d ), we see that 
(up to subsequences) there are a set Q C M d of unit measure such that ln n — > 1q in 
L 1 (R d ) and a function w G such th at w n n — >■ wsj in L 2 (R d ). Moreover, w > 

on M rf and {u> > 0} C By Proposition 2.4 and the inequality R{ w> oy < we 
conclude that the compactness (i) holds. 

(ii) limr._5.oo Q( r ) = ck G (0, 1). Let e > 0, then there exits r £ > 1/e such that for every 
R > r £ we have a — e < Q(R) < a. By the monotonicity of Q n (r) and the pointwise 
convergence to Q(r) we can find R e > r £ + 1/e and iV e such that 

a — 2e < Qn(R) < a + e, for every n> N £ and every r £ < R < R £ . 

By the definition of Qn the above inequality implies that there is a sequence x n G M d 
such that 

a — 3e < |Q n n < a + e for every n> N £ and every r £ < R < R £ . 

Defining 

A e n = n n n B re („ n ) and = O n \ Br £ (x„), 
we see that, thanks to the choice of R £ , 

d(A £ n ,B £ n ) >R £ -r £ > 1/e, 

(6.2) 

\\An\ - a\ + \\B n \ - (1 - a)\ < 8e and \Q n \ (A £ n U B £ n )\ < Ae . 

Up to substitute r £ and R £ with some r £ G {r £ ,r £ + y_-) and i? e G (i? E — y 7 ., _t £ ), we 
may suppose that 

H d - l (dB re {x n ) n n n ) +^- 1 (_£? i?£ (x n ) n o n ) < 2V_, 

and, as a consequence, 

P(/4 U ££) < P(n n ) + 2V_. (6.3) 

It remains to estimate the difference wq u — wu n , where U n := A^UB^. Let </> G C^°(M. d ) 
be positive with support in B 2 and equal to 1 on B\. For r > 0, consider the function 
4> r (x) = 4>{x/r). Defining 

U X n := <t> re /2{- ~ Xn)WQ n G H^{A £ n ), u\ : = (1 - (j) Re (- - X n )) W Un G Hq(B*), 

we have that 

\\wn n -u\- u 2 n \\ L 2 < 4_||u;n n ||oo, (6.4) 

where we noticed that we may choose r £ and R £ still satisfying all the previous as- 
sumptions and such that 

|(fi n n B 2Re (x n )) \ B rE/2 (x n )\ < 4e. 
Moreover, there is some universal constant C > such that 

I \Vu 1 n \ 2 dx+ [ \Vu 2 n \ 2 dx- I \Vw n J 2 dx < — < Ce. (6.5) 

JR d JR d JR d r e 
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where the last inequality follows by the choice of r e we made at the beginning of the 
proof. 

Since U n C Q n , we have that wu n is the orthogonal projection of WQ n on the space 
Ho(U n ) with respect to the Hq scalar product V Hence 



V(wn n -w Un )\ 2 dx < / \V(wn n -w* - u 2 n )\ 2 dx 



\Vw n J 2 dx-2 (ul + uDdx 



\S7u\\ 2 dx + 



\Vul\ 2 dx 



(6.6) 



J (wn n -ui- u 2 n ) dx + j (\Vui\ 2 + \Vu 2 n \ 2 - |V^J 2 ) dx < Ce. 



where in the first and second equality we have taken into account the equation satisfied 



by wn n while in the last inequality we have used (6.4) and (6.5). Sending e — > we 



see that (6.2) gives points (b) and (c) of the dichotomy case statement, (6.3) gives 



point (d), and (6.6) together with Proposition |2.5| gives (e). Since (a) is trivially true 
we obtain that in this case the dichotomy (ii) holds. 

(iii) lim^oo Q(r) = 0. In this case the first part of the statement of the vanishing case is 
clear, while for the second one we need some further considerations, similar to those 
in [6j Proposition 3.5], based on the Lieb's Lemma (see [28J). First notice that by a 
truncation argument, it is enough to prove the statement in the case when all £l n are 
bounded sets. Since \\Rn n ||£(L 2 (R d )) = -^i(^n) 1 a it is enough to prove that 



lim Ai(O n ) = +oo. 



(6.7) 



Let e > be fixed, R > be large enough and N S N be such that for all n > N and 
all i 6 R, we have \Q n n Br(x)\ < e. By the Lieb's Lemma, we have that there is 
some x G M. d such that 



Ai(0„ n B R (x)) < Ai(O n ) + X 1 (B R {x)). 

~ 2 Xi(Bi), for any r > 0, and that the ball minimizes Ai among 



Using that Xi(B r ) = 
all sets of given measure, we have 

MiBje-W - R~ 2 X 1 {B 1 ) < Ai(fi n ). 



(6.8) 



Since the left-hand side of (6.8) goes to infinity as e — > 0, we obtain (6.7). 



Let us assume now that (6.1) does not hold. In this case, clearly ln n — > in L 1 (M rf ), 



hence the same arguments of case (iii) imply that we are in the vanishing case. 
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